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Abstract. Chemical vapor deposition (CVD) is a widely used technique to 
grow solid materials with accurate control of layer thickness and composition. 
Under mass-transport-limited conditions, the surface of thin films thus produced 
grows in an unstable fashion, developing a typical motif that resembles the 
familiar surface of a cauliflower plant. Through experiments on CVD production 
of amorphous hydrogenated carbon films leading to cauliflower-like fronts, we 
provide a quantitative assessment of a continuum description of CVD interface 
growth. As a result, we identify non-locality, non-conservation, and randomness as 
the main general mechanisms controlling the formation of these ubiquitous shapes. 
We also show that the surfaces of actual cauliflower plants and combustion fronts 
obey the same scaling laws, proving the validity of the theory over seven orders of 
magnitude in length scales. Thus, a theoretical justification is provided, that had 
remained elusive thus far, for the remarkable similarity between the textures of 
surfaces found for systems that differ widely in physical nature and typical scales. 



1. Introduction 

Chemical vapor deposition (CVD) is a technique that is extensively used to grow films 
whose surfaces have controlled smoothness or composition [I| . Part of the generalized 
use of CVD to produce coatings or thin films is due to the fact that it can be used 
with almost all elements and with many compounds. Basically, CVD involves the film 
growth of a solid out from the aggregation of species that appear as the result of the 
reaction or decomposition of volatile precursors within a chamber. Chemical reactions 
occur in the vicinity of or at the surface of the solid, by-products being removed when 
needed. Here, we are interested in CVD as a technique that is capable to grow a 
surface under far from equilibrium conditions, yielding unstable rough surfaces that 
resemble the morphology of a familiar cauliflower plant [2]. We will refer to these 
surfaces as cauliflower-like fronts. 

Interestingly, not only growing thin films display this appealing cauliflower 
texture, but also many other natural patterns do. In general, these shapes, although 
easily recognizable, are not regular but present some self-similar or hierarchical 
structure within a characteristic sea of randomness. In this sense, cauliflower-like 
fronts rank among the most fascinating natural forms, in view of their simplicity and 
considering their diversity in origins and scales: they can be observed across length 
scales that range from tens of nanometers (surfaces of amorphous thin films [2]) up 
to hundreds of microns (turbulent combustion fronts [3]) and tens of centimeters (the 
familiar cauliflower plants). However, these morphologies being originated under non- 
equilibrium conditions, there is a lack of a general theoretical framework that can 
account for such a diversity and ubiquity. 

Another feature that makes CVD attractive as a benchmark to understand surface 
growth far from equilibrium is the possibility to formulate a physically motivated 
theory for interface dynamics, which incorporates the essential mechanisms that 
drive the process when this production technique is employed [4]. However, to our 
knowledge, a detailed comparison between such theory and the mentioned (fractal) 
cauliflower-like fronts is still lacking. This is remarkable in view of the wide interest 
that fractal geometry [5] has raised in the past, having been recognized to encode the 
morphological features of self-similar systems, namely those whose structure looks the 
same with independence of the scale of observation. Actually, many of the best-known 
fractals — as, for instance, computational models of biological morphogenesis [6] — are 
geometrical structures constructed deterministically by iteration of a simple initial 
motif. However, this qualitative knowledge is not entirely satisfactory because, as 
mentioned, cauliflower-like structures are not exactly regular but, rather, appear 
random to the eye. 

In this work we provide a detailed comparison between experimental surfaces of 
thin amorphous hydrogenated carbon films grown by CVD and predictions from a 
physical model derived from first principles. Excellent agreement is obtained both 
for (qualitative) morphological as well as for (quantitative) statistical analysis. This 
allows us to identify the main features of cauliflower-like fronts, as well as the essential 
general mechanisms that lead to their occurrence, thus accounting for their ubiquity 
in natural systems across several orders of magnitude. These conclusions are reached 
after further morphological analysis of actual cauliflower plants and combustion fronts 
for which typical scales are macroscopic, rather than submicrometric as in our CVD 
experiments. The interface evolution equation we consider is thus postulated as a 
universal description for non-local interface growth under appropriate conditions. 
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2. Model 



In Refs. 00 [9] a generic system of equations for CVD surface growth was presented. 
Those equations contain the main mechanisms involved in CVD growth: diffusion in 
the vapor phase, reaction, and attachment via surface kinetics, generalizing the classic 
description of the process (see [4] and references therein) to account for fluctuation 
effects in aggregation events and diffusive fluxes. Performing a standard linear and 
weakly nonlinear analysis, one arrives at a closed equation for the height of the film 
surface, h(r,t), at time t, where r is the position above a reference plane. Actually, 
the equation is more easily expressed for the space Fourier transform of the surface, 
hq(t) = ^[^(x, t)] and, using q = |q|, reads 

dtK = (Vq - Dd q 3 ) \ + j -F{(V^) 2 } q + 77 q . (1) 

Herdl V is the average velocity of the interface, D is the diffusion coefficient in 
the vapor phase, dg is the capillarity length of the surface, and the Gaussian white 
noise term rjq contains the information about the underlying microscopic fluctuations 
[3 [H [9] , having zero mean and correlations given by 

v7q(*toq'(0> =D n (2ir) 2 6(q + q')6(t-t'). (2) 

Despite the apparent simplicity of equation (JTJ) , emphasis must be done in its real space 
representation in order to stress its non-local character. Thus, the terms proportional 
to odd powers of q correspond to fractional Laplacians acting on the height field, 

(-VT- 1)/2 Mr)=c 2 , 2p _ lP v/ y ~ dr', p=l,2, (3) 

where PV denotes Cauchy principal value and C2,2p-i are appropriate numerical 
constants jTOJ [TTJ [J3] . Indeed, the Fourier representation of © is given by 

T[{ ^)(2p-l)/2 h]=q 2p-l h ^ p= l j2) (4) 

as occurring in Eq. (p}. Thus, in real space these terms couple height differences 
with algebraically decaying kernels. Hence, the value of the local growth velocity at 
a given surface point depends on the values of the height at all other surface points. 
Physically, this non-local coupling arises from the competition of the different parts 
of the system over the resources for growth [14] . In the case of a solid surface growing 
out of species that aggregate from a vapor phase, it is induced by the geometrical 
shadowing of prominent surface features that are more exposed to diffusive fluxes, 
over more shallow ones |4Ki 



3. Comparison with experiments 

In order to understand the physical implications of Eq. ([lj, we resort to numerical 
simulations that circumvent the analytical difficulties posed by its nonlinearity. 
Moreover, we have performed growth experiments in order to show that this 
equation indeed describes quantitatively actual cauliflower-like morphologies in CVD. 
Specifically, amorphous hydrogenated carbon (a-C:H) films were grown by electron 
cyclotron resonance chemical vapor deposition (ECR-CVD) on silicon substrates in a 
commercial ECR reactor (ASTeX, AX4500) in a two-zone vacuum chamber operating 

| See Table II for a glossary and summary of the main variables and parameters used in this work. 
§ Similarly, for combustion of premixed flames, the competition is for the available unburnt fuel |13| . 
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Experiments Theory 




Figure 1. Nanocauliflower surface growth. AFM top views (lxl /im 2 ) of CVD 
experiments at times t = 40 min (A), f = 2 hrs (B), and f = 6 hrs (C). Panels 
D-F show numerical results from Eq. {TJ for the same times. 
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Observable Value 

Growth velocity (V) 2.4 ± 0.12 x 10" s cm s" 1 (864 ± 43 nm/h) 

Mean substrate temperature (T) 343 K 

Mean free path (£ m f ) 0.45 ± 0.05 cm 

Partial pressure (methane) 3.75 ± 0.25 x 10 -3 mbar 

Diffusion coefficient (D) 0.33 ± 0.02 crn^ 1 



Table 1. Summary of parameters that can be measured or estimated from our 
CVD growth experiments. 



with a 2.45 GHz microwave source at 208 — 210 Watt input power. Gas mixtures 
of methane/argon (15 sccm/35 seem) were applied keeping the operating pressure 
at 1.1 x 10 -2 Torr. A dc bias of —50 V was applied to the silicon substrates. All 
samples were grown under these conditions and only the deposition time was varied. 
The film surface morphology was characterized ex-situ by Atomic Force Microscopy 
(AFM) with a Nanoscope Ilia equipment (Veeco) operating in tapping mode with 
silicon cantilevers. 

Top views of the surface morphology are shown in Fig. 1, panels A to C. As we 
see, a "globular" structure appears at short times with a characteristic length-scale, 
that grows in a disorderly fashion with further deposition. In our ECR-CVD growth 
system, the main growth species are ions and radicals. The latter can be distinguished 
into two main groups, i.e., CiH x and C2H X radicals. Within the first subgroup, 01112,3 
radicals have values of the sticking probability s to become permanently attached to 
the surface about s ~ 10~ 4 — 10~ 2 , whereas CiH and Ci have a sticking coefficient 
close to unity [15]. The C2H X radicals generally have a high sticking coefficient 
(s ~ 0.4 — 0.8) 16J. In fact, for a pure methane plasma, an overall sticking coefficient 
s = 0.65 ± 0.15 has been estimated [IT]. Moreover, when methane is diluted with 
argon, the impingement of argon ions generates dangling bonds at the surface leading 
to an effective increase of s for the different growth species. Thus, we can assume 
that the effective sticking coefficient is close to unity in our system, s ~ 1 [18 . This 
allows us to determine the surface kinetics regime at which experiments are operating. 
Assessment is done through comparison of the two velocity scales in the system: the 
mean surface velocity, V, and the mass transfer rate ko, that is related to the sticking 
probability s through 

kD=DL-] p ^-, (5) 

with L m { p being the mean free path of molecules in the vapor phase (see [19) and 
references therein). From data in Table 1 and Eq. ((5]), we can estimate ko — 0.75 cm 
s _1 , which is considerably higher than V — 2.4 x 10~ 8 cm s _1 . Hence, the system 
can be assumed to be in the fast kinetics regime for which Eq. ([IJ is expected to hold 

laisiEi- 

In spite of the previous assessment of parameter values, still many microscopic 
details of the experimental setup cannot be measured or even estimated from data, 
mainly due to the limited resolution of the experimental measurements, and also due 
to the coexistence between species both in the vapor phase and at the very same 
aggregate surface. To cite a few, the mean atomic volume of aggregating species at 
the surface, surface tension, or the capillarity length. Unfortunately, some of these 
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Name 


Meaning 


CVD 


Chemical Vapor Deposition 


ECR-CVD 


Electron Cyclotron Resonance CVD 


T 


Fourier transform operator 


V 


Average interface growth velocity 


D 


Diffusion coefficients of particles in the gas 


d 


Capillarity length 


lc 


Characteristic lengthscale in the linear regime 


D n 


Amplitude of the noise fluctuations 


C-2,3,4 


Stabilizing coefficients in the general model {e.g., C3 = d^D) 


h{x,t) 


Position (height) of the growing interface 


Kit) 


Fourier transformation of /i(x, t) 




Noise term accounting for fluctuations 


v*(t) 


Fourier transformation of ?y(x, t) 


PSD 


Power spectral density (also 5(q, t)) 


W(t) 


Global roughness or width if the interface (standard deviation of the height) 


a 


Roughness exponent 


z 


Dynamic exponent 


P 


Growth exponent (/3 = a/z) 



Table 2. Summary of the acronyms and main variables used in this work. 



are crucial in order to determine the quantitative values of the coefficients in Eq. ([I]). 

Hence, in order to proceed further we must extract additional parameter values 
from analysis of the morphologies in Fig. fA-C. First, we render Eq. ([T]) non- 
dimensional by fixing appropriate time, length, and height scales, namely, we perform 
the following change of variables: x — > x' = x/xq (so q — > q' = x$q, t — > t' = i/to, 
h — > h! = h/ho, so that Eq. ([T]) reads 

h dt<l = (h Vq > _ Dd ^( q r) K, + ^$H{v'ti?U> + (* - Vr v V- (6) 

Cq \Xq Xq J z,Xq 

By properly choosing x , t Q and ho, as 

[Ddv~ [Ddo~ 
x = h = ^—, t = ^—, (7) 

we can reduce the latter equation to (after dropping the primes for convenience) 

d t K = (q- q 3 ) + i^{(V^) 2 } q + (|^) ry q , (8) 

where ry q is a white noise term with zero mean and variance 1. Note how all the 
information is now contained in the prefactor of the noise term which is, after fixing the 
lateral size of the simulation domain (in our case L = 512, see the discussion below), 
the only remaining free parameter to be fitted. We have performed simulations for 
different values of D n until we have found the optimal value that provides the best 
agreement with the experiments (in our case, tQD n /h^x^ = 0.2f). 

Starting from a flat initial condition, at short times surface slopes are small so that 
the quadratic nonlinear term in Eq. (TTJ) is expected to be negligible. The system will 
thus evolve according to the linear terms, which has implications on the statistical 
properties of the morphology. For instance, the computed skewness of the height 
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distribution in Fig. 1A is negligible (note how the height distribution in Fig. 2A is 
almost symmetric), which is consistent with a small contribution of the nonlincarity 
|V/i| 2 , that is the only term breaking the up-down symmetry of the surface. Neglecting 
this term a characteristic length scale can be identified in Eq. (fT]) that, prior to non- 
dimensionalization, is given by 



Hence, the ratio between l c and the lateral system size L x in numerical simulations of 
Eq. (fT|) must agree with the ratio between the experimental value l c = 28 nm and the 
experimental AFM window, L c — 1 /im. Thus, we obtain approximately L x ~ 512 in 
our dimensionless units (we have rounded this value up to an exact power of 2 in order 
to optimize the numerical integration of the equation by means of a pseudo-spectral 
algorithm) . 

The lengthscale l c is the geometric average of the diffusion length (in the bulk), 
Id = D/V) and the capillarity length. Physically, this average arises from the 
competition between the scales explored by the diffusing particles in the gas and the 
lengthscales at which they can travel on the surface until the either aggregate or 
evaporate. In practice, l c can be interpreted as the typical size of the cauliflower-like 
structures than can be identified in the surface morphology at short times, see Fig. 



Numerical simulations of Eq. ([T]) using the same scheme as in [5] are shown in 
Fig. 1D-F for the same set of times as for the experimental images that appear in the 
same figure. The time evolution of the surface consists of an initial regime controlled 
by the Mullins-Sekerka linear instability [2Tj that leads to the appearance of a pattern 
(cusp) with characteristic length scale l c . In a process that is reminiscent of the 
stochastic Michelson-Sivashinsky equation that describes combustion fronts (see [2"U] 
and references therein), there is a competition between cusp coarsening/annihilation, 
and cusp formation induced by noise, leading to fully non-linear dynamics. As a result, 
for long enough times unstable growth is stabilized by the quadratic Kardar-Parisi- 
Zhang (KPZ) nonlinearity [20 , the surface morphology becoming disordered and 
rough, with height fluctuations that are scale free both in time and in space 7] . Note 
the strong resemblance between the theoretical and the experimental morphologies 
shown in Fig. 1. 

We have done a more quantitative comparison between Eq. (JIJ and experimental 
surfaces. In particular, we have determined the distribution of heights (Fig. 2), 
the height power spectral density (PSD, see Figs. 3 and 4), S(q,t), defined as 
S(q,t) — (hq(t)h—q(t)), where brackets denote average over noise realizations, and 
the global width or roughness, W(i) (Fig. 3), defined as the standard deviation of the 
interface height around its mean. 

The distribution of heights for short and long times is shown in Fig. 2 
corresponding to the morphologies shown in Figs. 1A,C,D,F. The distributions are 
very noisy but, overall, the shape of the curves is comparable for both experiments 
and theory. As mentioned above, for short times (t — 40 minutes) the system is in 
the linear regime and one can neglect the role of the non-linearity. As a result, all the 
terms preserve the symmetry h — > —h and the distributions are symmetric. On the 
other hand, for long times larger slopes develop as a result of the initial exponential 
growth and the non-linear term breaks that symmetry. This can be easily seen in Fig. 




(9) 



1A. 



2B. 
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Ah (|im) Ah (jxm) 



Figure 2. Normalized distribution of heights, P(Ah) for A) Figures 1A (black 
solid line) and ID (blue dashed line) corresponding to t = 40 minutes and B) 
Figures 1C (black solid line) and IF (blue dashed line) for t = 6 hours. 



As discussed in the introduction, cauliflower-like fronts are characterized by scale 
invariance, which can be quantified with the power spectral density (PSD). For the 
time scales of Fig. 1 , the PSD reflects the scale invariance associated with kinetic 
roughening (self-affine interfaces), and is expected to behave as S(q) ~ q-( 2 "+<i) |22) . 
Here, a is the so-called roughness exponent and d = 2 is the substrate dimension. 
Besides, the roughness grows as a power law of time W(t) ~ t a / z , where z is the so- 
called dynamic exponent that measures the speed at which height correlations spread 
laterally across the interface [22] . It is customary to define a third roughness exponent, 
/3, that is related to the previous ones through j3 = a/z. For each curve S(q,t) in 
Fig. 3, the small q behavior corresponds to an uncorrelated interface, the crossover 
to correlated spectra moving to smaller q (larger length scales) as time proceeds. In 
our case we obtain numerically a — 1.03 ± 0.06 and (3 — 0.93 ± 0.07 which are equal, 
within error bars, to the values a — (3 — z = 1 predicted by Renormalization Group 
(RG) calculations on Eq. ([1]) [2Ql [23] . In order to obtain these exponent values, we 
have computed the roughness using the relation 

W 2 (t) = J S(q,t)dq = J 2irqS(q,t)dq, (10) 

where S(q,t) is the radially averaged PSD (note that the argument here is q = |q|). 
Hereafter we will refer to this radially averaged function simply as 2D PSD. 

In Fig. 3A we show experimental time evolution of the surface roughness. As 
shown, the value obtained for the growth exponent f3 — 0.93 ± 0.07 is close to the 
theoretical prediction (8 = 1. A customary method to determine the roughness 
exponents is by means of the collapse of the PSDs at different times, by properly 
scaling S(q,t) and q as shown in Figs. 3B (experiments) and Fig. 3C (theory). 

The quantitative agreement that we obtain between the experimental and 
theoretical PSD functions is not limited to the values of the scaling exponents. 
Remarkably, it actually extends to the behavior of the full functions along the 
dynamics of the system, as can be seen in Fig. 4, where we compare theoretical 
and experimental PSDs for short and long times. Such type of agreement goes much 



8 




Figure 3. A) Experimental global width or roughness, W(t) (solid circles). The 
blue solid line has slope = 0.93 ± 0.07. B) Collapse of the radially averaged 
power spectral density for times from t = 40 minutes to 6 h (see legend) obtained 
for a = 1.03 and z = a//3, with f) as obtained from W(t). C) Same as panel B) 
but for the theoretical model. 




Figure 4. Scaling universality of cauliflower surface growth. Normalized PSD 
functions are shown to compare systems spanning several orders of magnitude in 
size. Relevant scales are indicated with arrows. The label "Thin Films" identifies 
the comparison between the experimental and theoretical (angular average of the) 
PSD for thin films grown by CVD. Circles: experiment t = 40 min. Solid red line: 
theory t = 40 min. Squares: experiment t = 6 h. Dashed blue line: theory 
t = 6 h. The orange dashed straight line is a guide to the eye with slope —4 
as predicted by Eq. |[T}. The solid black line under the label "Cauliflowers" has 
been obtained after averaging results obtained for 10 cauliflower slices (Brassica 
Oleracea, from two different specimens). The solid brown line under "Combustion 
fronts" corresponds to the PSD of an experimental combustion profile |24| . The 
green dashed lines have slopes —3. This value differs from the —4 obtained for 
the Thin Films case because we are computing the PSD of one-dimensional slices 
in the cases of the cauliflowers and combustion fronts. 



beyond what is usually expected on the context of universal properties [22 
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4. Universality of cauliflower- like fronts 

Thus far we have been using the term cauliflower in a loose way. In order to justify this 
usage it is convenient to identify the main universal phenomena and morphological 
properties that occur in the emergence of cauliflower-like structures in other contexts, 
and compare them with those assessed in the previous section. For instance, in the 
case of the familiar plants one can postulate (i) an interaction among the branches that 
sustain the external surface. This interplay would induce competitive growth among 
different plant features; moreover (ii) mass is non conserved and (iii) fluctuations 
are intrinsic to the biological underlying processes taking place both at the level 
of the cell metabolism and in the interaction with the environment, (iv) An extra 
stabilizing ingredient is necessary in order to guarantee the dynamical stability of the 
ensuing surface (whose specific form, as we argue below, does not change the statistical 
properties of that surface). In line with the occurrence of universality in the properties 
of rough surfaces evolving far from equilibrium [22) , for appropriate cases the argument 
can be reversed. Thus, when comparing two different systems, the same statistical 
properties of the surface are a manifestation of the same governing general principles, 
in spite of the fact that the detailed physical mechanisms controlling the dynamics of, 
e.g., aggregating species in CVD and plant cells, are quite different indeed. 

This property is a generalization of what happens in the proximity of a critical 
point within the framework of critical phenomena in equilibrium systems. Thus, while 
the microscopic details are different, the character of the interactions (in our case 
imposed by non-local competition or non-conservation) dictates the dynamics. For 
instance, in the case of the cauliflowers, non-locality is caused by branch competition. 
In the case of CVD, it stems from the fact that the diffusing particles access with 
a higher probability the most exposed parts of the surface. Finally, in the case of 
combustion, the parts of the front which lie behind the average have less access to 
oxygen and other combustion species. 

The basic ingredients expected for cauliflower-like surface growth, (i) through (iii) 
above, should reach a non-trivial balance resulting generically into a morphology that, 
albeit disordered, presents a self-similar, hierarchical structure. Moreover, one expects 
a typical characteristic length-scale to arise at the finest observation scale, due to the 
competition between stable and unstable growth mechanisms, as generically occurs in 
pattern forming systems [21] , 

Schematically, 

Variation of 
height 

Non-conservation 
(non-linearity) 

If we were to formulate a general interface equation that incorporates these 
mechanisms, actually arguments exist that can impose restrictions on its possible 
mathematical form. First, we assume that such an equation is weakly non-linear, 
in the sense that it is a polynomial in small powers of the height and its derivatives. 
This is a standard simplification in the study both of scale- invariant spatially extended 
systems [25], and of pattern forming systems [5T] in the long wavelength limit close 
to instability threshold, although see [53] and below. Next, self-similarity requires 
system statistics to remain unchanged under amplification of the lateral length scale 



Non-locality 
(competition) 



Stabilizing mechanism 
(short scale) 



Fluctuations 
(noise) 



(11) 
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of the sample by a factor b, while rescaling the surface height by the same factor, 
(h(br)) = (bh(r)), where (•) means in a statistical sense. For instance, in the case of 
cauliflower plants, one would expect the exposed surface to arise as an envelope for 
an underlying branched structure. This branching structure provides volume-filling 
mechanisms that guarantee efficient distribution of energy and nutrients. 

Moving further, the dominant linear term in the sought-for equation of motion 
can be inferred with large generality through dimensional analysis. Thus, assuming 
there is a single velocity scale, V, involved in the surface growth, the (linear) rate of 
amplification, u\, of a fluctuation can be related with the typical length-scale of the 
perturbation, A, as 



This expression is traditionally referred to as a dispersion relation, and is often 
written in terms of the wave number q = |q| = 2ir\~ 1 . Actually, Eq. (|12[) ensues 
for the celebrated Diffusion-Limited Aggregation model (DLA) that is the paradigm 
of fractal growth [27], and contains the signatures of unstable growth and non-local 
branch competition. In addition, further stabilizing mechanisms contribute to Eq. ()12j) 
as higher powers of q. Physical examples of non-local growth include solidification 
from a melt [28], flame fronts [29], stratified fluids [30], thin film evolution due to 
crystalline stress [31] , viscous fluid fingering [32] , growing biomorphs [33] , or geological 
structures [34], to cite a few. 

The final ingredient for the height equation of motion is non-linearity. The 
natural choice is the KPZ term (V/2)(Vh) 2 , that has been argued to be generically 
present in the continuum description of surfaces that grow irreversibly in the absence of 
conservation laws [35] and has been recently assessed to a high degree of accuracy in 
one-dimensional experiments [36] • Likewise, the simplest expression of fluctuations 
is through a random (uncorrelated) function of space and time like a Gaussian 
distributed white noise ~q(Y,t). 

Combining all these general ingredients together, we can write down the evolution 
equation for the local surface velocity, that in view of Eq. (| 12[) takes a particularly 
simple form when written for the Fourier modes of the height and noise fields, namely, 



where Cj are negative constants that depend on the specific stabilizing physical 
conditions. Note that Eq. (TTJ) simply corresponds to the particular case of Eq. (fTB"]) 
in which C2 = C4 = 0. An important result is that the same values of the scaling 
exponents a = z = 1 occur for any stabilizing linear mechanism of the form Cjq 3 
with j > 2, as indicated by RG analysis [3D]. Consequently, we are confident that 
this scaling behavior can be also identified in other systems for which the stabilizing 
term may have different non-zero contributions Cj. For instance, combustion fronts 
[29] and stratified fluids [30] correspond to C2 7^ and C3 = 0, and one again finds 
a = z = 1. 

In practical terms, all the terms in the main equation cooperate to produce those 
self-similar structures. Thus, if we replace the term Vq for, for instance, vq 2 then 
one would have obtained the celebrated Kuramoto-Sivashinsky equation. Or, if one 
suppresses the nonlinear term (^-J-[Wh) 2 ] then the numerical integration will explode 
as the surface roughness would increase exponentially without control. 



ui\W = dimensionless constant =>• uj q ~ Vq. 



(12) 




4 



(13) 



j=2 



11 



As an example, in Fig. 4 we also show the power spectral density of the profile 
of an expanding spherical flame (the experimental profile was taken from Ref. [24]). 
showing good agreement with the predictions of Eq. (JT3J) - Additionally, we have 
performed statistical measurements of the outermost surfaces of cauliflower plants. 
Specifically, we have computed the power spectral density of the interface profiles of 
several slices from specimens of Brassica Oleracea. The roughness exponent a can 
be determined by analyzing the power spectral density of these slices. In Fig. 4 we 
also show the averaged PSD over 10 different slices. The slope for small values of 
q leads to the value a = 1.02 ± 0.05 that, within experimental uncertainty [37] > is 
the same as the one provided by our continuum description. Overall, Fig. 4 proves 
the validity of the generalized Eq. (fl"3|) to describe cauliflower-like fronts across seven 
orders of magnitude in length scales, finally justifying the use of the term "cauliflower" 
as applied to morphologies that can differ quite strongly from the familiar plants. The 
difference between these and CVD films or combustion fronts is that, for the former, 
we have not been able to obtain a dynamical characterization of the plant morphology 
(we have only characterized it at a fixed, long time), as opposed to the latter in which 
the full dynamical equation Eq. (Ti~3|) has been derived from first principles and has 
been also experimentally validated, both in the context of CVD (this work) and for 
combustion fronts, see [38] . 

To better understand the significance of the unit values found for both critical 
exponents a and z, note that for kinetically rough, self-affine surfaces, W/L ~ L a ~ l 
for a sufficiently large lateral observation scale, L [22]. Precisely for a — 1 the system 
is not merely self-affine but becomes, rather, self-similar, its geometrical features 
remaining statistically invariant in the macroscopic limit L —> oo . Note moreover that 
in this case the average growth velocity for Eq. (I13[) is only due to the nonlinear term 
and becomes scale-independent precisely for a = 1, showing the self-consistency of 
our initial assumption on a single velocity scale. Also the fact that z = 1 reflects 
another peculiar fact about the fractality of the system: the system is self-similar 
also in time. Namely, correlations travel ballistically across the surface so that time 
behaves as space under rescaling. Hence, if we observe the system at two different 
times, we cannot distinguish the second one from an isotropic spatial zoom performed 
in the earlier one, as illustrated in Fig. [5] This corresponds to the intuition that by 
mere visual inspection it is hard to distinguish between the whole cauliflower plant 
and a piece of it, and between young and small florets. Therefore, in spite of the lack 
of dynamical information about cauliflower plant growth, the fact that self-similarity 
in time constraints the value of z to be unity, and the confirmation of a ~ 1 from Fig. 
3, both give us confidence to suggest that the growth of cauliflower plants obeys the 
same general principles as CVD growth. 

5. Conclusions 

To conclude, we have shown that Eq. (jTJ) provides an accurate description of 
unstable thin film growth by CVD, agreeing with experiments both qualitatively 
as well as quantitatively. To our knowledge, this is the first time in which such a 
quantitative agreement between theory and experiment is achieved, that goes beyond 
values of critical exponents, reaching the full dynamical behavior of observables such as 
the power spectral density. Actually, the moving boundary problem leading to Eq. (TTJ) 
has quite a generic form that is relevant to a number of processes in which transport is 
diffusive, such as solidification from a melt or electrochemical deposition (see references 
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Figure 5. Fractality in space and time reflects the inability to distinguish between 
image zooms and time shifts. This is the hallmark of cauliflower-like surface 
growth as described by Eq. 1131 . Top: Numerical simulation for ti = 2 h. Bottom 
right: Numerical simulation for = 4 h. Bottom left: Spatial amplification of 
the top panel by a factor t2/t\ = 2. The zoomed area is indicated in the top 
panel with a white dashed square. 

e.g. in [20]), so that similar quantitative descriptions of different growth systems by 
Eq. {T]) can be foreseen. Beyond that, we have also seen that a similar scaling behavior 
characterized by scaling exponents a = z = 1 can be moreover described by Eq. (|13p 
that applies to other systems that differ in the (linear) relaxation mechanisms. 

We would like to emphasize some important points concerning the implications of 
Eq. (|13[) . It is intriguing that the geometrical properties of cauliflower- like structures 
are at the boundary between disorder and fractality, between self-affinity and self- 
similarity. Thus, the values a = z — 1 of the scaling exponents induce an interface 
which is disordered at all scales, while allowing at the same time for the identification 
of a "typical" texture or motif. RG calculations and numerical simulations [3UJ H3] 
both indicate the robustness of these exponents values, suggesting the universality of 
Eq. (TT3|) as a description of a large class of non-equilibrium systems. Note, however, 
that interfaces developed under the same general physical principles as elucidated here, 
but for which the evolution equation is strongly, rather than weakly, nonlinear, may 
feature different morphological properties from the present cauliflower type. Examples 
are known in the dynamics of thin |39j and epitaxial films, and are reviewed in |26j . 

One of the reasons why fractals are so popular is the promise that, knowing 
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their generating rules, we can infer the character of the underlying physical or 
biological mechanisms. Hence whether the interactions are non-local vs local, 
non-conserved vs conserved, self-similar vs self-affine, . . . dictates the form of the 
mathematical equations. In contrast to "algorithmic" descriptions of fractals, the 
virtue of our continuum dynamical formulation is that it allows us to extract which 
are the most relevant mechanisms [40] whose interplay gives rise to these appealing 
structures, namely, non-locality, non-conservation, and noise. Among all the possible 
mathematical forms of non- locality, self-similarity enforces a = 1. This conclusion is 
expected to guide the inference of the relevant mechanisms at play in specific physical 
or biological systems where cauliflower-like structures are identified. Moreover, it 
is remarkable that such a simple equation as Eq. (II 3[) can be able to capture this 
non-trivial dynamics, to the extent that, by means of pseudo-spectral numerical 
integration, the system is capable of efficiently producing realistic patterns that 
resemble turbulent flame fronts or the texture of cauliflower plants. 

From a more general point of view, our theory also brings up the long-standing 
question as to why natural evolution favors self-similar structures. The so-called 
allometric scaling relations 41 explain (and predict) the branching structure of living 
bodies. The central idea behind these theories is that biological time scales are limited 
by the rates at which energy can be spread to the places where it is exchanged with the 
tissues. Thus, the space-filling structure 42 required to supply matter and energy to 
a living system can be accounted for. Focusing on more specific systems (cauliflower 
plants, etc.), albeit with a large degree of universality, our work suggests the self- 
similar features that the "canopy" atop such branched structures may have. 
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